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>. Abstract 

,_l ' The highly radiopure ~ 250 kg Nal(Tl) DAMA/LIBRA set-up is running at 

■^- , the Gran Sasso National Laboratory of the I. N.F.N. . In this paper the first result 

£ — ■ obtained by exploiting the model independent annual modulation signature for 

("N| Dark Matter (DM) particles is presented. It refers to an exposure of 0.53 tonxyr. 

^-4- , The collected DAMA/LIBRA data satisfy all the many peculiarities of the DM 

f****; ■ annual modulation signature. Neither systematic effects nor side reactions can 

QQ | account for the observed modulation amplitude and contemporaneously satisfy all 

the several requirements of this DM signature. Thus, the presence of Dark Matter 

|L| particles in the galactic halo is supported also by DAMA/LIBRA and, considering 

■ H ! the former DAMA/Nal and the present DAMA/LIBRA data all together (total 

j^ ■ exposure 0.82 tonxyr), the presence of Dark Matter particles in the galactic halo 

$H ' is supported at 8.2 a C.L.. 

Keywords: Scintillation detectors, elementary particle processes, Dark Matter 

PACS numbers: 29.40.Mc - Scintillation detectors; 95.30.Cq - Elementary particle 
processes; 95.35.+d - Dark matter (stellar, interstellar, galactic, and cosmological) . 

1 Introduction 

DAMA/LIBRA is part of the DAMA project, which is mainly based on the develop- 
ment and use of low background scintillators [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 
15]. 

In particular, the former DAMA/Nal and the present DAMA/LIBRA set-ups have 
the main aim to perform a direct detection of Dark Matter (DM) particles in the 



galactic halo through the model independent annual modulation signature (originally 
suggested in rcf. [16]). 

The former DAMA/Nal experiment [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14] has 
achieved many competitive results on rare processes and, in particular, has pointed 
out a model independent evidence for the presence of DM particles in the galactic halo 
with high C.L.. 

In 1996 - while running the DAMA/Nal set-up - DAMA proposed to INFN to 
develop and build a one ton set-up [17] to further investigate Dark Matter particles and 
other rare processes. Thus, a second generation R&D project was funded to develop 
new highly radiopure Nal(Tl) detectors towards the achievement of an intermediate 
step: the ~ 250 kg highly radiopure Nal(Tl) DAMA/LIBRA (Large sodium Iodide 
Bulk for RAre processes) set-up, which is now in data taking. 

The exploitation of the annual modulation DM signature with highly radiopure 
Nal(Tl) as target material can permit to answer - by direct detection and in a way 
largely independent on the nature of the candidate and on the astrophysical, nuclear 
and particle Physics assumptions - the main question: "Are there Dark Matter (DM) 
particles in the galactic halo?" 

In particular, the use of the highly radiopure DAMA/LIBRA (and, previously, 
DAMA/Nal) Nal(Tl) scintillators as target-detectors offers many specific advantages 
thanks e.g. to the intrinsic radiopurity, to the large sensitivity to many of the DM 
candidates, of the interactions and of astrophysical, nuclear and particle Physics sce- 
narios, to the granularity of the set-up, to the data taking up to the MeV scale (even 
though the optimization is made for the lowest energy region), to the full controll of 
the running conditions, etc.. 

Phcnomcnological properties of some basic interaction mechanisms induced by DM 
particles are discussed, for instance, inref. [16, 18, 19, 20, 4, 5, 6, 8, 9, 10, 11]. The DM 
annual modulation signature exploits the effect of the Earth revolution around the Sun 
on the number of events induced by DM particles in a suitable low background set-up 
placed deep underground. In particular, as a consequence of its annual revolution, 
the Earth should be crossed by a larger flux of DM particles around roughly June 2 nd 
(when its rotational velocity is summed to the one of the solar system with respect 
to the Galaxy) and by a smaller one around roughly December 2 nd (when the two 
velocities are subtracted) [4]. Thus, the contribution of the signal to the counting rate 
in the fc-th energy interval can be written as (see e.g. ref. [4, 5]): 

Sk = S ,k + S m . k cosuj{t - t ) , (1) 

where: i) Sq k is the constant part of the signal; ii) S m ,k is the modulation amplitude; 
hi) u) = ?£■ with period T; iv) to is the phase. 

This annual modulation signature is very distinctive since a seasonal effect induced 
by DM particles must simultaneously satisfy all the following requirements: 1) the 
rate must contain a component modulated according to a cosine function; 2) with 
one year period; 3) with a phase roughly around June 2 nd in case of usually adopted 
halo models (slight variations may occur in case of presence of non thermalized DM 
components in the halo); 4) this modulation must be present only in a wcll-dchned 
low energy range, where DM particles can induce signals; 5) it must be present only 
in those events where just a single detector, among all the available ones in the used 



set-up, actually "fires" (single-hit events), since the probability that DM particles 
experience multiple interactions is negligible; 6) the modulation amplitude in the region 
of maximal sensitivity has to be < 7% in case of usually adopted halo distributions, but 
it may be significantly larger in case of some particular scenarios such as e.g. those of 
ref. [19]. To mimic such a signature spurious effects or side reactions should be able not 
only to account for the observed modulation amplitude but also to contemporaneously 
satisfy all the requirements of the signature; none of these has been found or suggested 
by anyone over more than a decade (see e.g. ref. [4, 5], the references therein and 
later) . 

The corollary question: "Which are exactly the nature of the DM particle (s) de- 
tected by the annual modulation signature and the related astrophysical, nuclear and 
particle Physics scenarios?" requires subsequent model dependent corollary analyses 
as those performed e.g. in refs. [4, 5, 6, 7, 8, 9, 10, 11]. One should stress that it 
does not exist any approach to investigate the nature of the candidate in the direct 
and indirect DM searches which can offer these information independently on assumed 
astrophysical, nuclear and particle Physics scenarios. 

As regards complementary information from accelerators, and most noticeably from 
LHC, the existence of some of the possible candidates could be tested at some extent; 
this will be very useful. However, it is worth noting that interesting DM candidates 
and scenarios for them exist, which are beyond the reach of that class of experiments 
(but potentially capable of determining the annual modulation effect). 

The main goal of the DAMA/LIBRA experiment is to further study the presence 
of DM particles in the galactic halo pointed out by the former DAMA/Nal experiment 
[4, 5] exploiting the annual modulation signature and to get improved information 
on the corollary quests on the nature of the candidate particle(s) and on the related 
astrophysical, nuclear and particle physics models. Moreover, second order effects 
are planned to be investigated (see e.g. ref. [5, 7]), and dedicated data takings will 
also allow the study of many other rare processes (as e.g. already performed with 
DAMA/Nal [12, 13, 14]) thanks to the peculiarity of the experimental set-up. 

DAMA/LIBRA and the former DAMA/Nal are the only experiments effectively 
exploiting in all the aspects the DM annual modulation signature, and with highly 
radiopure Nal(Tl) detectors. Note that approaches based on many selections and 
handling procedures to "reject" the elect romagntic component of the counting rate 
cannot offer any signature for Dark Matter particles even under the assumption of an 
"ideal" electromagnetic component rejection, since e.g. the neutrons and the internal 
end-range a's induce signals indistinguishable from recoils (they are looking for) which 
cannot be estimated and subtracted in any reliable manner at the needed precision, 
and since part or all the signal can have electromagnetic nature. Moreover, in a safe 
investigation of the DM annual modulation signature those data handlings cannot be 
applied e.g. because of their always statistical nature which would affect the annual 
modulation analysis and restrict the sensitivity to many kinds of candidates (including 
also the WIMPs). On the other hand, as known, the exploitation of the DM annual 
modulation signature acts itself as an effective background rejection. 

In this paper some main features of the DAMA/LIBRA set-up [21] will be shortly 



summarized in sect. 2. The model independent experimental results obtained by 
DAMA/LIBRA (exposure of 0.53 tonxyr collected over 4 annual cycles) and the com- 
bined ones with DAMA/Nal (exposure of 0.29 tonxyr collected over 7 annual cycles) 
are presented (total exposure of 0.82 tonxyr) in sect. 3. In sect. 4 the quantitative 
investigation on possible systematic effects and side processes is discussed. The corol- 
lary model dependent analyses on the candidate particle(s) and astrophysical, nuclear 
and particle physics scenarios will be presented elsewhere in a dedicated publication; 
here in Appendix A just few arguments are mentioned for some illustrative purposes. 



2 The experimental set-up 

The DAMA/LIBRA set-up, its main features and radiopurity have been discussed in 
the devoted ref. [21]. Here we just shortly summarize some information. 

The installation of DAMA/LIBRA started in July 2002 after the dismounting of 
the former DAMA/Nal. The experimental site as well as many components of the 
installation itself have been implemented. All the procedures performed during the 
dismounting of DAMA/Nal and the installation of DAMA/LIBRA detectors have been 
carried out in high purity (HP) Nitrogen atmosphere. 

The sensitive part of DAMA/LIBRA is made of 25 highly radiopure Nal(Tl) crystal 
scintillators in a 5-rows 5-columns matrix. Each Nal(Tl) detector has 9.70 kg mass and 
a size of (10.2 x 10.2 x 25.4) cm 3 . The bare crystals are enveloped in Tetratec-tcflon 
foils and encapsulated in radiopure OFHC Cu housing; 10 cm long special quartz light 
guides act also as optical windows on the two end faces of the crystals and are coupled 
to two low background photomultipliers (PMT). The threshold of each one of the two 
PMTs on a detector is set at single photoelectron level; their coincidence provides 
the trigger of the detector. The software energy threshold has been cautiously taken 
at 2 keV electron equivalent (hereafter keV). The measured light response is 5.5-7.5 
photoelectrons/keV depending on the detector. The detectors are housed in a low ra- 
dioactivity sealed copper box installed in the center of a low-radioactivity Cu/Pb/Cd- 
foils/polyethylene/parafhn shield; moreover, about 1 m concrete (made from the Gran 
Sasso rock material) almost fully surrounds (mostly outside the barrack) this passive 
shield, acting as a further neutron moderator. The copper box is maintained in HP 
Nitrogen atmosphere in slightly overpressure with respect to the external environment; 
it is part of the 3-levels sealing system which excludes the detectors from environmen- 
tal air. The whole installation is air-conditioned and the temperature is continuously 
monitored and recorded; moreover, it is worth noting that the detectors have copper 
housings in direct contact with the multi-tons passive shield and its huge heat ca- 
pacity (w 10 6 cal/°C ) further assures a relevant stability of the detectors operating 
temperature (see also later). 

Following the same strategy as DAMA/Nal, on the top of the shield a glove-box 
(also continuously maintained in the HP Nitrogen atmosphere) is directly connected 
to the inner Cu box, housing the detectors, through Cu pipes. The pipes are filled 
with low radioactivity Cu bars (covered by 10 cm of low radioactive Cu and 15 cm of 
low radioactive Pb) which can be removed to allow the insertion of radioactive sources 
for calibrating the detectors in the same running condition, without any contact with 



external air. The glove-box is also equipped with a compensation chamber. 

A hardware/software system to monitor the running conditions is operative and 
self-controlled computer processes automatically control several parameters and man- 
age alarms. For the electronic chain, the data acquisition system and for all other 
details see ref. [21]. 

The DAMA/LIBRA set-up, as the former DAMA/Nal, allows the recording both 
of the single-hit events (those events where just one detector of many actually fires) 
and of the multiple-hit events (those events where more than one detector fire). 

The experiment take data up to the MeV scale despite the optimization is made 
for the lowest energy region. The linearity and the energy resolution of the detectors 
at low and high energy have been investigated using several sources as discussed in 
ref. [21]. In particular, as regards the low energy region, calibrations down to the 3.2 
keV X-ray have been carried out . During the production runs periodical calibrations 
(every ~ 10 days) are carried out with 241 Am sources, introduced in the proximity of 
the detectors by source holders inserted in the Cu pipes mentioned above; the latter 
one is also continuously maintained in the HP Nitrogen atmosphere. 

The energy threshold, the PMT gain, the electronic line stability are continuously 
verified and monitored during the data taking by the routine calibrations, by the 
position and energy resolution of internal lines [21] and by the study of the hardware 
rate behaviours with time. 

The main procedures of the DAMA data taking for the investigation of DM par- 
ticles annual modulation signature are: 1) the data taking of each annual cycle starts 
from autumn/winter (when cosu>(t — i ) — 0) towards summer (maximum expected); 
2) the routine calibrations with radioactive sources are performed about each 10 days 
(collecting typically ~ 10 4 — 10 5 events per keV), moreover regularly intrinsic calibra- 
tion are carried out, etc. [21]; 3) the on-line monitoring of all the running parameters 
is continuously carried out with automatic alarm to operator if any would go out of 
allowed range. 



3 The model-independent experimental results 

As mentioned, DAMA/LIBRA started the first operations in March 2003. However, 
in order to allow the decay of medium half- life isotopes [21], the data taking for the 
investigation of the annual modulation signature, reported here, has been started on 
September 9, 2003. Moreover, it is worth noting that one of the more external detectors 
has been put out of operation few months after installation because of a PMT break; 
since the disinstallation and reinstallation of this detector would require the opening 
of the set-up, the installation of the room for doing it in HP Nitrogen atmosphere and 
some time of stop of the experiment, we have delayed this. The related procedure are 
planned to occur in 2008 when also an upgrade of the electronics will occur. Therefore, 
the exposed mass in the four annual cycles, presented here, is 232.8 kg for a total 
exposure of 0.53 tonxyr. 

The only data treatment, which is performed on the raw data, is to eliminate 
obvious noise pulses (mainly PMT noise, Chcrcnkov light in the light guides and in the 
PMT windows, and afterglows) near the energy threshold in the single-hit events [21]; 



the number of such pulses sharply decreases when increasing the number of available 
photoelectrons. In particular, as mentioned, the DAMA/LIBRA detectors are seen 
by two PMTs working in coincidence and this already strongly reduces the noise near 
the energy threshold for the single-hit events (of interest for the detection of DM 
particles), while obviously noise is practically absent in the multiple-hit events, since 
the probability to have random coincidences is negligible (< 3 x 10~ 6 ). This rejection 
of the noise near energy threshold is based on the different time structures of the pulse 
profile of noise pulses (time decay of order of tens ns) and of the scintillation pulses 
(time decay of order of hundreds ns). The high number of photoelcctrons/keV assures 
a very good separation between the two populations, nevertheless stringent acceptance 
windows are used in order to assure full noise rejection near energy threshold; related 
acceptance window efficiencies are measured by devoted source calibrations. For a 
description of the used procedure and details see ref. [21]. 

Detailed information about the four annual cycles by DAMA/LIBRA. considered 
here, is given in Table 1. In these annual cycles about 4.4 x 10 7 events have also 



Table 1: DAMA/LIBRA annual cycles. There a — (cos 2 ui(t — to)) is the mean value of 
the squared cosine and f3 = (cosu)(t — to)) is the mean value of the cosine (the averages 
are taken over the live time of the data taking and to = 152.5 day, i.e. June 2 nd ); thus, 
a — f3 2 indicates the variance of the cosine (i.e. it is 0.5 for a full year of data taking). 
The exposed mass in these four annual cycles is 232.8 kg; see text. 



Period 




Exposure (kg x day) 


a-/3 2 


DAMA/LIBRA- 1 


Sept. 9, 2003- July 21, 2004 


51405 


0.562 


DAMA/LIBRA-2 


July 21, 2004- Oct. 28, 2005 


52597 


0.467 


DAMA/LIBRA-3 


Oct. 28, 2005 - July 18, 2006 


39445 


0.591 


DAMA/LIBRA-4 


July 19, 2006 - July 17, 2007 


49377 


0.541 


Total 




192824 
~ 0.53 tonxyr 


0.537 



been collected for energy calibrations and about 6.0 x 10 7 events for the evaluation of 
the acceptance windows efficiency for noise rejection near energy threshold [21]. The 
periodical calibrations and, in particular, those related with the acceptance windows 
efficiency mainly affect the duty cycle of the experiment; in the present data taking it 
is of the same order as the one of DAMA/Nal, despite the larger number of involved 
detectors, thanks also to the improvements in the electronics and in the data acquisition 
(DAQ) system. 

Figure 1 shows the cumulative low-energy distribution of the single-hit scintillation 
events (of interest for the DM particles since DM particle multiple interaction proba- 
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Figure 1: Cumulative low-energy distribution of the single-hit scintillation events (that 
is each detector has all the others as veto), as measured by the DAM A/LIBRA detec- 
tors in an exposure of 0.53 ton x yr. The energy threshold of the experiment is 2 keV 
and corrections for efficiencies are already applied. 



bility is negligible), as measured by the DAMA/LIBRA detectors in the 0.53 tonxyr 
exposure. 

In order to further investigate the presence of DM particles in the galactic halo, a 
model-independent investigation of the annual modulation signature has been carried 
out by exploiting the time behaviour of the residual rates of the single-hit events in the 
lowest energy regions of the DAMA/LIBRA data. These residual rates are calculated 
from the measured rate of the single-hit events (obviously corrections for the overall 
efficiency and for the acquisition dead time are already applied) after subtracting the 
constant part: < r^fc — flatjk >jk- Here r^k is the rate in the considered z-th time 
interval for the j-th detector in the fc-th energy bin, while flatjk is the rate of the 
j-th detector in the fc-th energy bin averaged over the cycles. The average is made 
on all the detectors (j index) and on all the 1 keV bins (k index) which constitute 
the considered energy interval. The weighted mean of the residuals must obviously be 
zero over one cycle. 

Figure 2 shows the time behaviour (over three energy intervals) of the model- 
independent residual rates for single-hit events collected by the new DAMA/LIBRA 
experiment over four annual cycles (0.53 tonxyr). Those measured over seven annual 
cycles by the former DAMA/Nal experiment [4, 5] (0.29 tonxyr) are shown as well; 
the cumulative exposure of the two experiments is 0.82 tonxyr. The advantage of the 
increased exposed mass and exposure is evident. 

In particular, the residual rates in the (2-4), (2 - 5) and (2 - 6) keV energy intervals 
are depicted in Fig. 2; the experimental points present the errors as vertical bars and 
the associated time bin width as horizontal bars. The superimposed curves represent 
the cosinusoidal functions behaviours Acosuj(t — to) with a period T = — — 1 y r 
and with a phase to = 152.5 day (June 2 nd ) and modulation amplitudes, A, obtained 
by best fit over the whole data (DAMA/Nal & DAMA/LIBRA). The dashed vertical 
lines correspond to the maximum of the signal (June 2 nd ), while the dotted vertical 
lines correspond to the minimum. We note that, for simplicity, in Fig. 2 the same time 
binning already considered e.g. in ref. [4, 5] has been used. The result of this approach 
is similar by choosing other time binnings, as it is also evident from the analysis on 



2-4 keV 



ft 

I 

p. 



0.1 

0.08 







il (0 ( 29 tonx; 

ass *= 873 kg 






<- 








M\ 


! (tatget id 

I li 


) 


| 


(target mass = 232:8 kg) 

A A Js ^ 






jVjIi 


f 


, : i 






I 







500 1000 1500 2000 2500 



2-5 keV 



4000 4500 

Time (day) 



<-DA]yiA/LIBJ{A (0.53jtonxyr>* 




1000 1500 2000 2500 



2-6 keV 



3000 5500 4000 4500 

Time (day) 



I 

a 



0.1 
0.08 



-DAMA/Nal (0 ( 29 tonxyr) 
(target mass i= 873 kg) 




<-DAMA/LIBJJA (0.53jtonxyr>* 



1500 2000 2500 



3000 3500 4000 4500 

Time (day) 



Figure 2: Model-indcpcndcnt residual rate of the single-hit scintillation events, mea- 
sured by the new DAMA/LIBRA experiment in the (2-4), (2 - 5) and (2 - 6) keV 
energy intervals as a function of the time. The residuals measured by DAMA/Nal and 
already published in ref. [4, 5] are also shown. The zero of the time scale is January 
I s * of the first year of data taking of the former DAMA/Nal experiment. The exper- 
imental points present the errors as vertical bars and the associated time bin width 
as horizontal bars. The superimposed curves represent the cosinusoidal functions be- 
haviours Acosu>(t — to) with a period T = — = 1 yr, with a phase to — 152.5 day 
(June 2 nd ) and with modulation amplitudes, A, equal to the central values obtained by 
best fit over the whole data, that is: (0.0215 ± 0.0026) cpd/kg/keV, (0.0176 ± 0.0020) 
cpd/kg/kcV and (0.0129 ±0.0016) cpd/kg/keV for the (2 - 4) keV, for the (2 - 5) keV 
and for the (2 - 6) keV energy intervals, respectively. See text. The dashed vertical 
lines correspond to the maximum of the signal (June 2 nd ), while the dotted vertical 
lines correspond to the minimum. The total exposure is 0.82 tonxyr. 



modulation amplitudes given in the following. 

Table 2 summarizes the results obtained by fitting with the function A cos to (t — to)'. 
i) only the DAMA/Nal data [4, 5]; ii) only the DAMA/LIBRA data; iii) all the data 
together. A clear modulation is present in all the energy intervals and the periods and 
phases agree with those expected in the case of a DM particle induced effect. 



Tabic 2: Results obtained from the time behaviours of the residual rates of the single- 
hit scintillation events, collected by DAMA/Nal, by DAMA/LIBRA and by the two 
experiments all together in the (2-4), (2 - 5) and (2 - 6) keV energy intervals. The 
data have been fitted with the function: A cos ui(t — to). The last column shows the 
C.L. obtained from the fitted modulation amplitudes. See comments in the text. 





A (cpd/kg/keV) 


T = *!■ (yr) 


t (day) 


C.L. 


DAMA/Nal 
(2-4) keV 
(2-5) keV 
(2-6) keV 


0.0252 ±0.0050 
0.0215 ±0.0039 
0.0200 ±0.0032 


1.01 ±0.02 
1.01 ±0.02 
1.00 ±0.01 


125 ±30 
140 ± 30 
140 ± 22 


5.0(7 

5.5(7 
6.3(7 


DAMA/LIBRA 

(2-4) keV 
(2-5) keV 
(2-6) keV 


0.0213 ±0.0032 
0.0165 ±0.0024 
0.0107 ±0.0019 


0.997 ±0.002 
0.998 ± 0.002 
0.998 ±0.003 


139 ±10 
143 ±9 
144 ±11 


6.7(7 

6.9(7 

5.6(7 


DAMA/NaI+ DAMA/LIBRA 
(2-4) keV 
(2-5) keV 
(2-6) keV 


0.0223 ±0.0027 
0.0178 ±0.0020 
0.0131 ±0.0016 


0.996 ±0.002 
0.998 ± 0.002 
0.998 ±0.003 


138 ±7 
145 ±7 
144 ±8 


8.3(7 

8.9(7 
8.2<7 



It is worthwhile remarking how the larger exposed mass per annual cycle has im- 
proved the fit; for example, the x 2 / 'd.o.f. of the best fit of the (2 - 6) keV single-hit 
residual rate from DAMA/Nal plus DAMA/LIBRA given in Table 2 is 53.2/64. The 
period and phase substantially agree with T =1 yr and to = 152.5 day. 

We note that the difference in the (2 - 6) keV modulation amplitude between 
DAMA/Nal and DAMA/LIBRA depends mainly on the rate in the (5 - 6) keV en- 
ergy bin. In particular, the modulation amplitudes for the (2 - 6) keV energy in- 
terval, obtained when fixing exactly the period at 1 yr and the phase at 152.5 days, 
are (0.019 ± 0.003) cpd/kg/keV and (0.011 ± 0.002) cpd/kg/keV for DAMA/Nal and 
DAMA/LIBRA, respectively; thus, their difference is about ~ 2ct which correspond 
to a modest, but non negligible probability. This is further supported by the analyses 
of the modulation amplitudes of each single year of DAMA/Nal and DAMA/LIBRA 
experiments, as reported in Fig. 3. There the central values obtained by best fit over 
the whole data set (see Fig. 2) are also depicted. The x 2 test (\ 2 — 4.9, 3.3 and 8.0 
over 10 d.o.f. for the three energy intervals, respectively) and the run test (lower tail 
probabilities of 74%, 61% and 11% for the three energy intervals, respectively) accept 
the hypothesis at 90% C.L. that the modulation amplitudes are normally fluctuat- 
ing around their best fit values. Thus, the cumulative result from DAMA/Nal and 
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Figure 3: Modulation amplitudes of each single year of DAMA/Nal and 
DAMA/LIBRA experiments in the (2-4), (2 - 5) and (2 - 6) keV energy inter- 
vals. The same time scale as in Fig. 2 is adopted. The solid horizontal lines shows the 
central values obtained by best fit over the whole data set (see Fig. 2). The \ 2 test 
and the run test accept the hypothesis at 90% C.L. that the modulation amplitudes 
are normally fluctuating around the best fit values. See text. 

In conclusion, the DAMA/LIBRA data are in substantial agreement with those of 
DAMA/Nal and the cumulative analysis favours the presence of a modulated cosine- 
like behaviour at 8.2 a C.L. (see Table 2). Moreover, the x 2 test on the residual rates 
disfavours the hypothesis of unmodulated behaviour (A = 0) giving probabilities of 
1.3 x 1CT 4 {x 2 /d.o.f. = 117.7/67), 1.9 x 1CT 4 (x 2 /d.o.f. = 116.1/67) and 1.8 x 1(T 4 
(x 2 / d.o.f. = 116.4/67) for the three energy intervals, respectively. 

The same data of Fig. 2 have also been investigated by a Fourier analysis (per- 
formed according to ref. [22] including also the treatment of the experimental errors 
and of the time binning); in particular, Fig. 4 shows the result for the DAMA/LIBRA 
and for the cumulative exposure; the one for DAMA/Nal alone has been given in ref. 
[4, 5]. Here a clear peak for a period of 1 year is evident in the lowest energy interval 
(2-6) keV. 

In the same Fig. 4-left there is also shown the absence of modulation in the 
region (6 - 14) keV just above the region where the modulation is present in the 
DAMA/LIBRA data; in Fig. A-right the same is shown for the cumulative DAMA/Nal 
and DAMA/LIBRA data (for the DAMA/Nal data alone sec ref. [4, 5]). 
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Figure 4: Power spectrum of the measured single-hit residuals for the (2-6) keV 
(solid lines) and (6-14) keV (dotted lines) energy intervals calculated according to 
ref. [22], including also the treatment of the experimental errors and of the time bin- 
ning. The data refer to: left -just to the DAMA/LIBRA data; right - to the cumulative 
DAMA/Nal and DAMA/LIBRA data; the case of DAMA/Nal has been given in ref. 
[4, 5]. As it can be seen, the principal mode present in the (2-6) keV energy interval 
corresponds to a frequency of 2.705 x 1CT 3 d _1 and 2.737 x 10~ 3 d _1 , respectively 
(vertical lines); that is, they correspond to a period of ~ 1 year. A similar peak is not 
present in the (6-14) keV energy interval just above. 
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Figure 5: Experimental residuals in the (2 - 6) keV region and those in the (6 - 14) 
keV region just above for the DAMA/LIBRA data considered as collected in a single 
annual cycle. The experimental points present the errors as vertical bars and the 
associated time bin width as horizontal bars. The initial time is taken at August 7 th . 
The clear modulation is present in the lowest energy interval, while it is absent just 
above. See text. 



Fig. 5 compares the residuals in the (2 - 6) keV region and those in the (6 - 14) 
keV region just above for the DAMA/LIBRA data considered as collected in a single 
annual cycle. A clear modulation is present in the lowest energy interval, while it is 
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absent just above. In fact, the best fitted modulation amplitude in the (6 - 14) keV 
energy region is well compatible with zero: (0.0009 ± 0.0011) cpd/kg/keV. 

Finally, Fig. 6 shows - for various energy intervals - the experimental single-hit 
residual rates, as collected in a single annual cycle, for the total exposure of 0.82 
tonxyr (i.e. DAMA/Nal & DAMA/LIBRA). 
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Figure 6: Experimental single-hit residual rates, as collected in a single annual cycle, 
for the total exposure of 0.82 tonxyr (i.e. DAMA/Nal plus DAMA/LIBRA) in the (2 
- 4) keV (a), (2 - 5) keV (b), (2 - 6) keV (c) and (6 - 14) keV (d) energy intervals. 
The experimental points present the errors as vertical bars and the associated time 
bin width as horizontal bars. The initial time is taken at August 7 th . Fitting the data 
with a cosinusoidal function when fixing exactly the period at 1 yr and the phase at 
152.5 days, the following amplitudes are obtained: a) (0.0204±0.0026) cpd/kg/keV; b) 
(0.0166 ±0.0020) cpd/kg/keV; c) (0.0125 ± 0.0016) cpd/kg/keV; d) (0.0004 ±0.0010) 
cpd/kg/keV. Thus, a clear modulation is present in the lowest energy regions, while 
it is absent just above. 

To verify absence of annual modulation in other energy regions and, thus, to also 
verify the absence of any significant background modulation, the energy distribution 
measured during the data taking periods in energy regions not of interest for DM 
detection have been investigated. In fact, the background in the lowest energy region is 
essentially due to "Compton" electrons, X-rays and/or Auger electrons, muon induced 
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events, etc., which are strictly correlated with the events in the higher energy part 
of the spectrum. Thus, if a modulation detected in the lowest energy region would 
be due to a modulation of the background (rather than to a signal), an equal or 
larger (sometimes much larger) modulation in the higher energy regions should be 
present. For this purpose, also in the present case we have investigated the rate 
integrated above 90 keV, R90, as a function of the time. In Fig. 7 the distribution of the 
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Figure 7: Distribution of the percentage variations of R90 with respect to the mean 
values for all the detectors in the DAM A/LIB RA-1 to -4 annual cycles (histogram); 
the superimposed curve is a gaussian fit. See text and sect. 4. 

percentage variations of R90 with respect to the mean values for all the detectors during 
the D AM A/LIBRA- 1 to -4 annual cycles is given. It shows a cumulative gaussian 
behaviour with u ~ 1%, well accounted by the statistical spread expected from the 
used sampling time. This result excludes any significant background variation (see 
also later). 

Moreover, fitting the time behaviour of R90 with phase and period as for DM par- 
ticles, a modulation amplitude compatible with zero is found in each running period: 
-(0.05±0.19) cpd/kg, -(0.12±0.19) cpd/kg, -(0.13 ±0.18) cpd/kg, and (0.15±0.17) 
cpd/kg for D AM A/LIBRA- 1 to -4 annual cycles, respectively. This excludes the pres- 
ence of a background modulation in the whole energy spectrum at a level much lower 
than the effect found in the lowest energy region for the single-hit events; in fact, 
otherwise - considering the R90 mean values - a modulation amplitude of order of 
tens cpd/kg, that ~ 100cr far away from the measured value, would be present. Other 
arguments are also given in sect. 4. 

Finally, a further relevant investigation has been done by applying the same hard- 
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ware and software procedures, used to acquire and to analyse the single-hit residual 
rate, to the multiple-hits one. In fact, since the probability that a DM particle interacts 
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Figure 8: Experimental residual rates over the four DAM A/LIBRA annual cycles for 
single-hit events (open circles) - class of events to which DM events belong - and 
for multiple-hits events (filled triangles) - class of events to which DM events do not 
belong - in the (2-4), (2-5) and (2-6) keV energy intervals. They have been obtained 
by considering for each class of events the data as collected in a single annual cycle and 
by using in both cases the same identical hardware and the same identical software 
procedures. The initial time of the scale is taken on August 7 th . The experimental 
points present the errors as vertical bars and the associated time bin width as horizontal 
bars. See text. 
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in more than one detector is negligible, a DM signal can be present just in the single- 
hit residual rate. Thus, this allows the test of the background behaviour in the same 
energy interval of the observed positive effect. We remind that similar investigations 
have already been performed for the last two annual cycles of the DAMA/Nal experi- 
ment [5]. Thus, in Fig. 8 the residual rates of the single-hit events measured over the 
four DAMA/LIBRA annual cycles are reported, as collected in a single annual cycle, 
together with the residual rates of the multiple-hits events, in the considered energy 
intervals. 

While, as already observed, a clear modulation is present in the single-hit events, 
the fitted modulation amplitudes for the multiple-hits residual rate are well compati- 
ble with zero: -(0.0004 ± 0.0008) cpd/kg/keV, -(0.0005 ± 0.0007) cpd/kg/kcV, and 
-(0.0004 ± 0.0006) cpd/kg/keV in the energy regions (2-4), (2 - 5) and (2 - 6) 
keV, respectively. Thus, again evidence of annual modulation with proper features as 
required by the DM annual modulation signature is present in the single-hit residuals 
(events class to which the DM particle induced events belong), while it is absent in 
the multiple-hits residual rate (event class to which only background events belong). 
Since the same identical hardware and the same identical software procedures have 
been used to analyse the two classes of events, the obtained result offers an additional 
strong support for the presence of a DM particle component in the galactic halo further 
excluding any side effect either from hardware or from software procedures or from 
background. 

The annual modulation present at low energy can also be shown by depicting the 
S m .k values as a function of the energy; the S m ^k is the modulation amplitude of the 
modulated part of the signal (see cq. 1) obtained by maximum likelihood method over 
the data considering T —I yr and to = 152.5 day. For such purpose the likelihood 
function of the single-hit experimental data in the k— th energy bin is defined as: 

N iik 

Lk = n ije -^^ (2) 

where Nijk is the number of events collected in the i-th time interval (hereafter 1 day), 
by the j-th detector and in the fc-th energy bin. N^k follows a Poissonian distribution 
with expectation value /i^fe = [bjk + £o,fc + S m ,k ■ cosio(ti — to)] MjAtiAEejk- The 
bjfc are the background contributions, Mj is the mass of the j— th detector, Ati is 
the detector running time during the i-th time interval, AE is the chosen energy 
bin, ejk is the overall efficiency. The usual procedure is to minimize the function 
jjk = — 2ln(Lu) — const for each energy bin; the free parameters of the fit are the 
(bjk + S'o.fc) contributions and the S m ± parameter. 

In Fig. 9 the S mt k are reported for the seven annual cycles of DAMA/Nal and for 
the four annual cycles of DAMA/LIBRA in each considered energy bin (here AE = 0.5 
keV). It can be inferred that positive signal is present in the (2-6) keV energy interval, 
while S m .k values compatible with zero are present just above. In fact, the S m .k values 
in the (6-20) keV energy interval have random fluctuations around zero with x 2 equal 
to 24.4 for 28 degrees of freedom. All this confirms the previous analyses. 
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Figure 9: Energy distribution of the S mi k variable for the total exposure (0.82 tonxyr, 
DAMA/Nal & DAM A/LIBRA). See text. A clear modulation is present in the lowest 
energy region, while S m ,k values compatible with zero are present just above. In fact, 
the Sm t k values in the (6-20) keV energy interval have random fluctuations around 
zero with \ 2 equal to 24.4 for 28 degrees of freedom. See also Appendix A. 



The method also allows the extraction of the S m (hereafter the index k is omitted) 
values for each detector, for each annual cycle as well as for each considered energy 
bin. The S m are expected to follow a normal distribution in absence of any systematic 

effects. Therefore, in order to show if they are statistically well distributed in all the 

s — ( s ) 
crystals, in all the annual cycles and in the energy bins, the variable x = — - — L -^ is 

considered. Here, a are the errors associated to S m and (S m ) are the mean values of the 
S m averaged over the detectors and the annual cycles for each considered energy bin 
(in the following AE = 0.25 keV). Similar investigations have already been performed 
also for DAMA/Nal [4, 5]. 

Figure 10 shows the distributions of the variable x for the DAM A/LIBRA data in 
the (2-6) keV energy interval plotted for each detector separately (i.e. the entries of 
each histogram are the 64 x values, evaluated for the 16 energy bins in the consid- 
ered (2-6) keV energy interval and for the 4 DAMA/LIBRA annual cycles). These 
distributions allow one to conclude that the observed annual modulation effect is well 
distributed in all the detectors and annual cycles. In fact, the standard deviations 
of the x variable for the DAMA/LIBRA detectors range from 0.80 to 1.16 (see also 
Fig. 11-bottom). Defining \ 2 — ^ x ~~ where the sum is extended over all the 64 x 
values - x 2 /d.o.f. values ranging from 0.7 to 1.28 (see Fig. 11 -top) are obtained. The 
corresponding upper tail probabilities range from about 97% to 6%. Therefore, the 
observed annual modulation effect is well distributed in all the detectors at 95% C.L.. 
The x 2 /d-o.f. values of the DAMA/LIBRA detectors show a distribution around their 
expectation value (see Fig. 11-top). The twenty- four points follow a \ 2 distribution 
with 64 degrees of freedom; in fact, when compared with the expectation in Fig. 11 - 
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Figure 10: Distributions (histograms) of the variable — ^ — K —^- (where a is the error 
associated to the S m and (S m ) are the mean values of the modulation amplitudes 
averaged over the detectors and the annual cycles for each considered energy bin) . The 
superimposed curves are gaussian fits. Each panel refers to a single DAMA/LIBRA 
detector in the (2 - 6) keV energy interval for the four DAMA/LIBRA annual cycles. 

top-right, a \ 2 /d.o.f. = 8.1/7 is obtained. The mean value of the twenty-four x 2 /d.o.f. 
values is 1.072, slightly larger than expected. Although this can be still ascribed to 
statistical fluctuations (see before), let us ascribe it to a possible systematics. In this 
case, one would have an additional error of < 5 x 10~ 4 cpd/kg/keV, if quadratically 
combined, or < 7 x 10~ 5 cpd/kg/keV, if linearly combined, to the modulation am- 
plitude measured in the (2 - 6) keV energy interval. This possible additional error - 
< 4.7% or < 0.7%, respectively, of the DAMA/LIBRA modulation amplitude - can be 
considered as an upper limit of possible systematic effects (see also later the dedicated 
section). 

The above arguments demonstrate that the modulation amplitudes are statistically 
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Figure 11: Top left: x 2 /d.o.f. values of S m distributions around their mean value for 
each DAMA/LIBRA detector in the (2-6) keV energy interval for the four annual 
cycles. The line at \ 2 /d.o.f. = 1.31 corresponds to an upper tail probability of 5%. 
All the x 2 /d.o.f . values are below this line and, thus, at 95% C.L. the observed annual 
modulation effect is well distributed in all the detectors. Top right: distribution of 
the twenty-four points in the top left panel (solid histogram) compared with the x 2 
distribution with 64 degrees of freedom; the comparison gives a x 2 1 d.o.f . = 8.1/7. See 
text. Bottom: standard deviations of the x variable for the DAMA/LIBRA detectors; 
see also Fig. 10. See text. 



well distributed in all the crystals, in all the data taking periods and in the considered 
energy bins. 

Let us, finally, release the assumption of a phase to = 152.5 day in the maximum 
likelihood procedure to evaluate the modulation amplitudes, as performed above; that 
is, let us alternatively write the expectation values, /iyfc, in eq. (2) as: 

Hijk = [b 3 k + So,k + S m ,k cos u)(ti - t ) + Z mM sinw(ii - t )} M^t^Ee^ (3) 
= [bjk + S , k + Y m . k cos w(U - t*)} MjAtiAEejk . (4) 

Obviously, for signals induced by DM particles one would expect: i) Z m .^ ~ (because 
of the orthogonality between the cosine and the sine functions); ii) S m ,k — Ym,k] hi) 
t* ~ to — 152.5 day. In fact, these conditions hold for most of the dark halo models; 
however, it is worth noting that slight differences can be expected in case of possible 
contributions from non-thcrmalized DM components, such as e.g. the SagDEG stream 
[7] and the caustics [23]. The analysis has been performed considering the data of 
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the seven annual cycles of DAMA/Nal and the four annual cycles of DAMA/LIBRA 
all together. Fig. 12 -left shows the 2a contours in the plane (S m ,Z m ) for the (2- 
6) keV and (6-14) keV energy intervals and Fig. 12-right shows, instead, those in 
the plane (Y m ,t*). The best fit values for the (2-6) keV energy interval are (la 
errors): S m = (0.0122 ± 0.0016) cpd/kg/keV; Z m = -(0.0019 ± 0.0017) cpd/kg/keV; 
Y m = (0.0123 ± 0.0016) cpd/kg/keV; t* = (144.0 ± 7.5) day; while for the (6-14) keV 
energy interval are: S m = (0.0005 ± 0.0010) cpd/kg/keV; Z m = (0.0011 ± 0.0012) 
cpd/kg/keV; Y m = (0.0012 ± 0.0011) cpd/kg/keV and t* obviously not determined 
(see Fig. 12). These results confirm those achieved above by analysing the residuals. 
In particular, a modulation amplitude is present in the lower energy intervals and 
the period and the phase agree with those expected for DM induced signals. Finally, 
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Figure 12: 2a contours in the plane (S m , Z m ) (left) and in the plane (Y m , t*) (right) for 
the (2-6) keV and (6-14) keV energy intervals. The contours have been obtained by 
the maximum likelihood method, considering the seven annual cycles of DAMA/Nal 
and the four annual cycles of DAMA/LIBRA all together. 

forcing to zero the contribution of the cosine function in eq. (3), the Z m ^ values 
as function of the energy have also been determined by using the same procedure. 
Obviously, such values are expected to be zero in case of presence of a DM particles' 
signal with t Q = 152.5 day. From the data we get: -(0.0046 ± 0.0041) cpd/kg/keV, 
-(0.0036±0.0042) cpd/kg/keV, (0.0036±0.0034) cpd/kg/keV, and -(0.0015±0.0032) 
cpd/kg/keV, for (2-3), (3-4), (4-5), (5-6) energy intervals, respectively. Moreover, the 
X 2 test applied to the data in the (2-14) keV energy region (x 2 equal to 12.4 for 12 
degrees of freedom) supports the hypothesis that the Z rn ^ values are simply fluctuating 
around zero. 

In conclusion, an annual modulation satysfying all the requirements of the DM 
annual modulation signature has been observed also in DAMA/LIBRA with high C.L.; 
it credits cumulatively with DAMA/Nal an evidence for the presence of DM particles 
in the galactic halo at 8.2 a C.L.. 



19 



4 Further investigation on possible systematic ef- 
fects and side reactions 

Also the data of the first four annual cycles of DAMA/LIBRA as previously those 
of DAMA/Nal fullfill the requirements of the DM annual modulation signature and 
prcliminar investigation on absence of any significant systematics or side reaction effect 
is already present in some parts of the previous section; however, here the argument 
will be addressed in more details. 

Sometimes naive statements are put forwards as the fact that in nature several 
phenomena may show annual periodicity. It is worth noting that the point is whether 
they might mimic the annual modulation signature in DAMA/LIBRA, i.e. whether 
they might be not only quantitatively able to account for the observed modulation 
amplitude but also able to contemporaneously satisfy all the requirements of the DM 
annual modulation signature. This was deeply investigated in the former DAMA/Nal 
experiment (see e.g. ref. [4, 5] and refs. therein; no one has been either found or 
suggested by anyone so far) and will be further addressed in the following for the 
present DAMA/LIBRA data. 

Firstly, in order to continuously monitor the running conditions, several pieces of 
information are acquired with the production data and quantitatively analysed; note 
that information on technical aspects of DAMA/LIBRA has been discussed in ref. 
[21], where the sources of possible residual radioactivity have also been analysed. 

In particular, all the time behaviours of the running parameters, acquired with the 
production data, have been investigated. Table 3 shows the modulation amplitudes 



Table 3: Modulation amplitudes (1 a error) obtained by fitting the time behaviours of 
the main running parameters including a possible annual modulation with phase and 
period as for DM particles. These running parameters, acquired with the production 
data, are: i) the operating temperature of the detectors; ii) the HP Nitrogen flux in the 
inner Cu box housing the detectors; iii) the pressure of the HP Nitrogen atmosphere 
of the inner Cu box housing the detectors; iv) the environmental Radon in the inner 
part of the barrack from which the detectors are however excluded (see text and ref. 
[21] for details); v) the hardware rate above single photoelectron threshold. All the 
measured amplitudes are compatible with zero. 





DAMA/LIBRA-1 


DAMA/LIBRA-2 


DAMA/LIBRA-3 


DAMA/LIBRA-4 


Temperature 


-(0.0001 ±0.0061)°C 


(0.0026 ±0.0086)°C 


(0.001 ±0.015)°C 


(0.0004 ±0.0047)°C 


Flux 


(0.13 ±0.22) 1/h 


(0.10 ±0.25) 1/h 


-(0.07 ±0.18) 1/h 


-(0.05 ±0.24) 1/h 


Pressure 


(15±30)10~ 3 mbar 


-(13±25)10~ 3 mbar 


(22±27)10 -3 mbar 


(1.8±7.4)10~ 3 mbar 


Radon 


- (0.029 ± 0.029) Bq/m 3 


-(0.030 ±0.027) Bq/m 3 


(0.015 ± 0.029) Bq/m 3 


-(0.052 ±0.039) Bq/m 3 


Hardware rate 


-(0.20±0.18)10- 2 Hz 


(0.09±0.17)10- 2 Hz 


-(0.03±0.20)10- 2 Hz 


(0.15±0.15)Kr 2 Hz 
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obtained for each annual cycle when fitting the time behaviours of the values of the 
main parameters including a cosine modulation with the same phase and period as for 
DM particles. As can be seen, all the measured amplitudes are well compatible with 
zero. 

Let us now enter in more details. 



4.1 The temperature 

Since temperature at sea level varies along the year, sometimes it has been naively 
suggested that it can mimic the observed effect. 

It is worth noting that the full experiment is placed underground and works in 
an air-conditioned environment; moreover, the detectors have Cu housing in direct 
contact with the multitons metallic passive shield whose huge heat capacity definitively 
assures a relevant stability of the detectors' operating temperature [21]. Nevertheless 
the operating temperature is read out by a probe and stored with the production data, 
in order to offer the possibility of quantitative investigations (see also above). 

Specific information on the D AM A/LIBRA- 1 to 4 annual cycles can be derived from 
Fig. 13; no evidence of any operating temperature modulation has been observed, as 
also quantitatively reported in Table 3. 
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(T-<T>)/<T> 

Figure 13: Distribution of the relative variations of the operating temperature mea- 
sured during the D AM A/LIBRA- 1 to -4 annual cycles (histogram); the superimposed 
curve is a gaussian fit. The standard deviation is 0.4%. 

However, to properly evaluate the real effect of possible variations of the detectors' 
operating temperature on the light output, we consider the distribution of the root 
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mean square temperature variations within periods with the same calibration factors 
(typically ~ 10 days); this is given in Fig. 14 cumulatively for the four-year data sets. 
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Figure 14: Distribution of the root mean square (r.m.s.) detectors' operating temper- 
ature variations within periods with the same calibration factors (typically ~ 10 days) 
during the DAM A/LIB RA-1 to -4 annual cycles. The mean value is 0.04 °C. 

Considering the obtained mean value of the root mean square detectors' operating 
temperature variation: ~ 0.04 °C, and the known value of the slope of the light output 
around its value: < -0.2%/°C, the relative light output variation is < 10~ 4 , which 
would correspond to a modulation amplitude < 10~ 4 cpd/kg/keV (that is < 0.5% of 
the observed modulation amplitude). 

Thus, there is no significant effect from possible temperature variation in the 
DAMA/LIBRA installation; on the other hand, in that case some of the specific re- 
quirements of the DM annual modulation signature (such as e.g. the 4 th and the 5 th ) 
would fail, while they are instead satisfied by the DAMA/LIBRA production data. 

In conclusion, all the arguments given above quantitatively exclude any role of 
possible effects on the observed rate modulation directly correlated with temperature. 

For the sake of completeness, we comment that sizeable temperature variations in 
principle might also induce variations in the electronic noise, in the Radon release from 
the rocks and in some environmental background; these specific topics will be further 
addressed in the following. 

4.2 The background 

As quantitatively discussed in previous section, in order to verify the absence of any 
significant background modulation, the energy distribution measured during the data 
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taking periods in energy regions not of interest for DM detection has been investigated. 
The presence of background (of whatever nature) modulation is already excluded by 
the results - given in previous section - on the measured rate integrated above 90 kcV, 
R90, as a function of the time; the latter one not only does not show any modulation, 
but allows one to exclude the presence of a background modulation in the whole 
energy spectrum at a level some orders of magnitude lower than the annual modulation 
observed in the single-hit events in the (2 - 6) keV energy region. 

A further relevant support is given by the result of the analysis of the multiple-hits 
events (see above) which independently proofs that there is no modulation at all in 
the background event in the same energy region where the single-hit events present 
an annual modulation satifying all the requirements of the DM annual modulation 
signature. 

The results given above obviously already account for whatever kind of background 
including that possibly induced by neutrons, by Radon and by side reactions. In the 
following we will focus our attention on the latter ones to offer independent cau- 
tious analyses to estimate their possible contribution, as done in ref. [3, 4] for the 
DAMA/Nal data. 

4.2.1 ... more on neutrons 

As regards the thermal neutrons, the neutron capture reactions 23 Na(n,7) 24 Na and 
23 Na(n,7) 24m Na (cross sections to thermal neutrons equal to 0.10 and 0.43 barn, re- 
spectively [24]) have been investigated. The 24 Na isotope is a /^-emitter (end point 
equal to 1.391 MeV) with two prompt associated 7's (2.754 and 1.369 MeV), while 
the 24m Na isotope decays 100% of the times in 24 Na by internal transition with a 7 of 
0.472 MeV. Thus, the possible presence of 24 Na isotopes in the Nal(Tl) crystals gives 
information about the thermal neutron flux surviving the DAMA/LIBRA shielding 
and impinging on the DAMA/LIBRA detectors; hence, as reported in ref. [21], this 
has been investigated with high sensitivity by looking for triple coincidences induced 
by a (3 in one detector and by the two 7's in two adjacent ones. An upper limit on the 
thermal neutron flux surviving the multicomponent DAMA/LIBRA shield has been 
derived as [21]: < 1.2 x 10~ 7 cm _2 s _1 (90% C.L.) 1 . The corresponding capture rate 
is: < 0.022 captures/day/kg. Assuming cautiously a 10% modulation (of whatever 
origin) of the thermal neutrons flux, the corresponding modulation amplitude in the 
lowest energy region has been calculated by MonteCarlo program to be < 0.8 x 10~ 6 
cpd/kg/kcV (that is < 0.01% of the observed modulation amplitude). In addition, 
a similar contribution cannot anyhow mimic the annual modulation signature since 
possible modulation of thermal neutron captures would induce e.g. variations in all 
the energy spectrum, that is it would fail some of the six requirements of the annual 
modulation signature (such as e.g. the 4 th and the 5 th ). 

A similar analysis has also been carried out for the fast neutrons case, as already 
done for DAMA/Nal [3, 4]. In particular, very safely, the effect of the about 1 m 



x We remind that the thermal neutron flux has been measured in the LNGS to be 1.08 X 10~ 6 
neutrons cm -2 s _1 [25] and two consistent upper limits on the thermal neutron flux have been 
obtained with the DAMA/Nal considering the same capture reactions and using different approaches 
[2, 13]. 
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concrete (made from the Gran Sasso rock material) which almost fully surrounds 
(outside the barrack) the DAMA/LIBRA passive shield acting as a further neutron 
moderator - is not cautiously included here. Thus, from the fast neutron flux measured 
at the Gran Sasso underground laboratory, ~ 10~ 7 fast neutrons cm~ 2 s _1 [26] the 
differential counting rate above 2 keV has been estimated by MonteCarlo code to be 
~ 10~ 3 cpd/kg/kcV. Therefore, assuming cautiously a 10% modulation (of whatever 
origin) of the fast neutron flux, the corresponding modulation amplitude is < 10~ 4 
cpd/kg/keV. 

Note that the use of other measurements of fast neutron flux at LNGS [25, 27] 
does not change the given conclusions. Moreover, an independent evaluation of the 
fast neutron flux impinging on the DAMA/LIBRA detectors has been performed by 
using the DAMA/LIBRA production data through the study of the inelastic reaction 
23 -/Va(n,n') 23 A r a*(2076 keV), which produces two j's in coincidence (1636 keV and 
440 keV). An upper limit - limited by the sensitivity of the method - has been found: 
< 2.2 x 10~ 7 fast neutrons cm~ 2 s^ 1 (90% C.L.), well compatible with the measured 
values in the laboratory in rcf. [25, 26, 27] and estimated in rcf. [28]. This further 
excludes any presence of a fast neutron flux in DAMA/LIBRA significantly larger than 
measured in ref. [25, 26, 27]. 

It is worth noting that a possible neutron flux modulation as claimed in rcf. 
[29] at LNGS (~ 5%) and in ref. [30] at the shallow deep Baradello mine or sug- 
gested by phenomenological approaches [28] 2 cannot quantitatively contribute to the 
DAMA/LIBRA observed effect, even if the neutron flux at LNGS would be assumed 
two orders of magnitude larger than measured. 

Finally, a possible modulation in the fast neutron flux would induce variation in 
all the energy spectrum and in the multiple-hits events at low energy, that is some of 
the six requirements mentioned above would fail. 

4.2.2 ... more on Radon 

The DAMA/LIBRA detectors are excluded from the air of the underground laboratory 
by a 3-level sealing system [21]; in fact, this air contains traces of the radioactive Radon 
gas ( 222 Rn - T1/2 = 3.82 days - and of 220 Rn - Ti/ 2 = 55 s - isotopes, which belong 
to the 238 U and 232 Th chains, respectively), whose daughters attach themselves to 
surfaces by various processes. In particular: i) the walls, the floor and the top of the 
inner part of the installation are insulated by Supronyl (permeability: 2 x 10~ n cm 2 /s 
[31]) and a large flux of HP Nitrogen is released in the closed space of this inner part 
of the barrack housing the set-up. An Oxygen level alarm informs the operator before 
entering it, when necessary; ii) the whole passive shield is sealed in a plexiglas box 
and maintained continuously in HP Nitrogen atmosphere in slight overpressure with 
respect to the environment as well as the upper glove box for calibrating the detectors; 
iii) the detectors are housed in an inner sealed Cu box also maintained continuously in 
HP Nitrogen atmosphere in slight overpressure with respect to the environment; the Cu 
box can enter in contact only with the upper glove box - during calibrations - which 



2 In particular, it has been suggested that a possible origin of the variability of the neutron flux 
might be ascribed to possible variations of the water content of the environment. 
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is also continuously maintained in HP Nitrogen atmosphere in slightly overpressure 
with respect to the external environment. 

Notwithstanding the above considerations, the Radon in the installation outside 
the plexiglas box, containing the passive shield, is continuously monitored; it is at level 
of sensitivity of the used Radonmeter. The time behaviours for the DAMA/LIBRA-1, 
-2, -3, and -4 annual cycles are shown in Fig. 15. As quantitatively reported in Table 
3, no modulation of Radon is present in the environment of the set-up; moreover, the 
detectors are further isolated by the other two levels of sealing [21]. 
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Figure 15: Time behaviours of the Radon in the inner part of the barrack (from which 
- in addition - the detectors are further isolated by other two levels of sealing [21]) 
during the DAMA/LIBRA-1 to -4 annual cycles, respectively. The measured values 
are at the level of sensitivity of the used radonmeter. The time scale has the origin at 
Jan. 1st, 2000. 

In Fig. 16 the distributions of the relative variations of the HP Nitrogen flux in 
the inner Cu box housing the detectors and of its pressure as measured during the 
DAMA/LIBRA-1 to -4 annual cycles are shown (the typical flux mean value for each 
annual cycle is of order of ~ 320 1/h and the typical overpressure mean value is of 
order of 3.1 mbar). 

We have also investigated possible Radon trace in the HP Nitrogen atmosphere 
inside the Cu box, housing the detectors, by searching for the double coincidences of 
the gamma-rays (609 and 1120 keV) from 214 Bi Radon daughter, obtaining an upper 
limit on the possible Radon concentration in the Cu box HP Nitrogen atmosphere: < 
5.8 x 10~ 2 Bq/m 3 (90% C.L.); thus, roughly < 2.5 x 10~ 5 cpd/kg/keV can be expected 
from this source in the lowest energy bins and single-hit events in DAMA/LIBRA as 
of interest for DM detection. The improvement - with respect to DAMA/Nal - of the 
limit of the expected contribution to the single-hit events counting rate at low energy 
is due to the enlarged matrix of detectors and to the better filling of the Cu box, 
housing the detectors, e.g. thanks to the new Cu shaped shield for PMTs and light 
guides [21]. 
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Figure 16: Distributions of the HP Nitrogen flux in the inner Cu box housing the 
detectors and of its pressure as measured during the D AM A/LIBRA- 1 to -4 annual 
cycles (histograms); the superimposed curves are gaussian fits. For clarity the HP 
Nitrogen flux has been given in terms of relative variations. 



This shows that even an hypothetical, e.g. 10%, modulation of possible Radon 
in the HP Nitrogen atmosphere of the Cu box, housing the detectors, would corre- 
spond to a modulation amplitude < 2.5 x 10~ 6 cpd/kg/keV (< 0.01% of the observed 
modulation amplitude). 

Moreover, it is worth noting that, while the possible presence of a sizeable quantity 
of Radon nearby a detector would forbid the investigation of the annual modula- 
tion signature (since every Radon variation would induce both the variation in the 
whole energy distribution and the continuous pollution of the exposed surfaces by 
the non- volatile daughters), it cannot mimic the DM annual modulation signature in 
experiments such as the former DAMA/Nal and the present DAMA/LIBRA which 
record the whole energy distribution; in fact, possible presence of Radon variation can 
easily be identified in this case and some of the six requirements of the DM annual 
modulation signature would fail. 

In conclusion, no significant effect is possible from the Radon. 



4.2.3 



more on side processes 



Finally, possible side reactions have also been carefully searched for. The only process 
which has been found as an hypothetical possibility is the muon flux modulation 
reported by the MACRO experiment [32]. In fact, MACRO has observed that the 
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muon flux shows a nearly sinusoidal time behaviour with one year period and maximum 
in the summer with amplitude of ~ 2 %; this muon flux modulation is correlated with 
the temperature of the Earth's atmosphere. A simple calculation to estimate the 
modulation amplitude expected from this process in the DAMA/LIBRA set-up can 
follow the analysis introduced in ref. [3, 4] and recalled in the following. The muon flux 
($u) and the yield of neutrons produced by muons measured at the underground Gran 
Sasso National Laboratory (Y) are: <J> M ~ 20 muons m _2 d _1 [32] and Y ~ (1— 7)xl0~ 4 
neutrons per muon per g/cm 2 [33] , respectively. Thus, the fast neutron rate produced 
by muons is given by: R n = <&^-Y -M e ff, where M e ff is the effective mass where muon 
interactions can give rise to events detected in the DAMA set-up. Consequently, the 
annual modulation amplitude in the lowest energy region induced in DAMA/LIBRA 
by a muon flux modulation as measured by MACRO [32] can be estimated according 
to: S$ = R n -g-e- f AE ■ f sing i e ■ 2%/(M set _ up ■ AE), where g is a geometrical factor, 
e is the detection efficiency for elastic scattering interactions, /ae is the acceptance 
of the considered energy window (E > 2 kcV), f single is the single-hit efficiency and 
2% is the MACRO measured effect. Since M set - up ~ 250 kg and AE ~ 4 keV, 
assuming the very cautious values g ~ e ~ /ae — /single — 0.5 and M e ff — 15 t, one 
obtains: Sm < (0.4 — 3) x 10~ 5 cpd/kg/keV. We stress that - in addition - the latter 
value has been overestimated of orders of magnitude both because of the extremely 
cautious values assumed in the calculation and, as mentioned, of the omission of the 
effect of the ~ 1 m concrete neutron moderation. Finally, we remark that not only 
the modulation of the muon flux observed by MACRO would give rise in our set-up 
to a quantitatively negligible effect, but - in addition - some of the six requirements 
necessary to mimic the annual modulation signature (such as e.g. the 4 th and the 5 th ) 
would fail. Therefore, it can be safely ignored. 

Just for the sake of completeness, we remind that the contribution of solar neutri- 
nos, whose flux is also expected to be modulated, is many orders of magnitude lower 
than the modulation amplitude measured by DAMA/LIBRA [14]. 

4.3 The noise 

Despite the good noise identification near energy threshold and the stringent noise 
rejection procedure which is used [21], the role of a possible noise tail in the data after 
the noise rejection procedure has been quantitatively investigated. 

In particular, the hardware rate of each detector above a single photoclectron, 
Rhj (j identifies the detector), has been considered. Indeed, this hardware rate is 
significantly determined by the noise. 

For the proposed purpose the variable: Rh = ^j{Rhj — {Rhj}), can be built; in the 
present case (Rhj) < 0.2 Hz. The time behaviours of Rh during the DAMA/LIBRA-1 
to -4 annual cycles are shown in Fig. 17. 

As can be seen in Fig. 18, the cumulative distribution of Rh for the DAMA/LIBRA- 
1 to -4 annual cycles shows a gaussian behaviour with a = 0.3%, value well in agreement 
with the one expected on the basis of simple statistical arguments. 

Moreover, by fitting the time behaviour of R# in the four data taking periods 
- including a modulation term as that for DM particles - a modulation amplitude 
compatible with zero: (0.03 ± 0.09) x 10 -2 Hz, is obtained. From this value the 
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Figure 17: Time behaviours of the variable Rh = T,j(Rhj — (Rhj)), where Rhj 
is the hardware rate of each detector above single photoelectron threshold (that is 
including the noise), j identifies the detector and (Rhj) is the mean value of Rhj in 
the corresponding annual cycle. 
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Figure 18: Distribution of R# during the DAMA/LIBRA-1 to -4 annual cycles (his- 
togram); the superimposed curve is a gaussian fit. See text. 
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upper limit at 90% C.L. on the modulation amplitude can be derived: < 1.8 x 1CP 3 
Hz. Since the typical noise contribution to the hardware rate of each detector is 
~ 0.10 Hz, the upper limit on the noise relative modulation amplitude is given by: 

1 ' S> 2 1 hHz IIZ - 7 - 2 x 1Cr4 ( 90% CX ')- Therefore, even in the worst hypothetical case of 
a 10% contamination of the residual noise - after rejection - in the counting rate, the 
noise contribution to the modulation amplitude in the lowest energy bins would be < 
7.2 x 10 -5 of the total counting rate. This means that an hypothetical noise modulation 
could account at maximum for absolute amplitudes less than 10~ 4 cpd/kg/keV. 

In conclusion, there is no evidence for any role of an hypothetical tail of residual 
noise after rejection. 

4.4 The calibration factor 

The performed calibrations have been discussed in ref. [21]; in particular, in long term 
running conditions periodical calibrations are performed every ~ 10 days with 241 Am 
source. 

Although it is highly unlikely that a variation of the calibration factor (proportion- 
ality factor between the area of the recorded pulse and the energy), tdcal, could play 
any role, a quantitative investigation on that point has been carried out. 

For this purpose, we define the percentage variation of each energy scale fac- 
tor (tdcal) with respect to the value measured in the previous calibration: e t dcai = 
" l tdcai ca k ~ 1 (here tdcalk is the value of the calibration factor in the fc-th calibra- 
tion). The distribution of e t dcai for all the detectors during the D AM A/LIBRA- 1 to 
-4 annual cycles is given in Fig. 19 -Left. This distribution shows a gaussian behaviour 
with a ~ 0.5%. Since the results of the routine calibrations are properly taken into 
account in the data analysis, such a result allows us to conclude that the energy cali- 
bration factor for each detector is known with an uncertainty <C 1% during the data 
taking periods. 

Moreover, the distribution of the percentage variations ((.he) of the high en- 
ergy scale factor with respect to the mean values for all the detectors and for the 
DAMA/LIBRA-1 to -4 annual cycles is reported in Fig. 19-right. Also this distribu- 
tion shows a gaussian behaviour with a ~ 0.6%. 

As discussed also in ref. [3, 4], the possible variation of the calibration factor for 
each detector during the data taking would give rise to an additional energy spread 
(o'cai) besides the detector energy resolution (a res ). The total energy spread can be, 
therefore, written as: a — \J a 2 es + <r 2 al ~ a res -[l + ^-(^ QSj -) 2 ]-, clearly the contribution 

due to the calibration factor variation is negligible since i • ( aca ,„ ) 2 < 7.5 x 10~ 4 „ n f,, 

(where the adimensional ratio „Jr v accounts for the energy dependence of this limit 
value). This order of magnitude is confirmed by a MonteCarlo calculation, which 
credits - as already reported in ref. [3, 4] - a maximum value of the effect of similar 
variations of tdcal on the modulation amplitude equal to 1 — 2 x 10~ 4 cpd/kg/keV. 
Thus, also the unlikely idea that the calibration factor could play a role can be safely 
ruled out. 
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Figure 19: Le/t: Distribution of the percentage variations (etdcai) of each energy scale 
factor (tdcal) with respect to the value measured in the previous calibration for the 
DAMA/LIBRA-1 to -4 annual cycles (histogram); the superimposed curve is a gaus- 
sian fit. The standard deviation is 0.5%. Right: Distribution of the percentage vari- 
ations (che) of the high energy scale factor with respect to the mean values for the 
DAMA/LIBRA-1 to -4 annual cycles (histogram); the superimposed curve is a gaus- 
sian fit. The standard deviation is 0.6%. 



4.5 The efficiencies 



The behaviour of the used overall efficiencies during the whole data taking periods has 
even been investigated. Their possible time variation depends essentially on the sta- 
bility of the efficiencies related to the adopted acceptance windows; they are regularly 
measured by dedicated calibrations [21]. In particular, Fig. 20 shows the percentage 
variations of the efficiency values in the (2-8) keV energy interval considering 2 keV 
bins. They show a gaussian distribution with a = 0.3% for DAMA/LIBRA-1 to -4, 
cumulatively. Moreover, we have verified that the time behaviour of these percentage 
variations does not show any modulation with period and phase expected for a possible 
DM signal. In Table 4 the modulation amplitudes of the efficiencies in each energy 
bin between 2 and 10 keV are reported, showing that they are all consistent with zero. 
In particular, modulation amplitudes - considering the four DAMA/LIBRA annual 
cycles all together - equal to (0.1 ± 0.4) x 10~ 3 and -(0.4 ± 0.4) x 10~ 3 are found 
for the (2-4) keV and (4-6) keV energy bins, respectively; both consistent with zero. 
Thus, also the unlikely idea of a possible role played by the efficiency values is ruled 
out. 
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Figure 20: Distribution of the percentage variations of the overall efficiency values 
with the respect to their mean values during the DAM A/LIB RA-1 to -4 annual cycles 
(histogram); the superimposed curve is a gaussian fit. See text. 



Table 4: Modulation amplitudes obtained by fitting the time behaviour of the efficien- 
cies including a cosine modulation with phase and period as for DM particles for the 
DAMA/LIBRA-1 to -4 annual cycles. 



Energy 



Amplitudes (xlO 



DAMA/LIBRA-1 



DAMA/LIBRA-2 



DAMA/LIBRA-3 



DAMA/LIBRA-4 



2-4 keV 
4-6 keV 
6-8 keV 
8-10 keV 



(0.3 ±0.6) 

(0.0 ±0.6) 

-(0.3 ±0.6) 

-(0.5 ±0.5) 



(0.1 ±0.6) 
-(0.7 ±0.6) 
-(1.0 ±0.7) 
-(0.5 ±0.5) 



-(0.4 ±1.1) 
-(0.3 ±1.0) 
-(0.2 ±0.8) 
-(0.2 ±0.6) 



-(0.4 ±1.0) 

-(0.7 ±1.0) 

-(1.0±0.8) 

(0.7 ±0.6) 



4.6 Conclusions on possible systematics effects 

No modulation has been found in any possible source of systematics or side reactions; 
thus, upper limits (90% C.L.) on the possible contributions to the DAMA/LIBRA 
measured modulation amplitude are summarized in Table 5. In particular, they cannot 
account for the measured modulation both because quantitatively not relevant and 
unable to mimic the observed effect. 
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Tabic 5: Summary of the results obtained by investigating possible sources of system- 
atics or of side reactions in the data of the DAMA/LIBRA-1 to -4 annual cycles. None 
able to give a modulation amplitude different from zero has been found; thus cautious 
upper limits (90% C.L.) on the possible contributions to the measured modulation 
amplitude have been calculated and arc shown here. 



Source 


Main comment 
(see also rcf. [21]) 


Cautious upper limit 
(90%C.L.) 


Radon 


Sealed Cu Box in 

HP Nitrogen atmosphere, 

3-level of sealing 


< 2.5 x ICT 6 cpd/kg/keV 


Temperature 


Air conditioning 
+ huge heat capacity 


< ICT 4 cpd/kg/keV 


Noise 


Efficient rejection 


< ICT 4 cpd/kg/keV 


Energy scale 


Routine 
+ intrinsic calibrations 


< I - 2 x 1(T 4 cpd/kg/keV 


Efficiencies 


Regularly measured 


< ICT 4 cpd/kg/keV 


Background 


No modulation above 6 keV; 

no modulation in the (2 - 6) keV 

multiple-hit events; 

this limit includes all possible 

sources of background 


< ICT 4 cpd/kg/keV 


Side reactions 


From muon flux variation 
measured by MACRO 


< 3 x ICT 5 cpd/kg/keV 


In addition: no effect can mimic the signature 



5 Conclusions 

The model independent results achieved by the second generation DAMA/LIBRA 
set-up in operation at the Gran Sasso National Laboratory confirms evidence of Dark 
Matter particles in the galactic halo with high confidence level; a cumulative C.L. of 8.2 
a is reached when considering the data of the former DAMA/Nal experiment and the 
present ones of DAMA/LIBRA all together. In particular, deep quantitative analyses 
exclude any effect either from systematics or from side processes (temperature, noise, 
hardware or software procedures, background of whatever nature including also radon, 
neutrons and cosmic rays). We note that no experiment exists whose result can be 
directly compared with those presented here. 

The wide sensitivity of the used target-detector material to many of the possible 
DM candidates and of the possible astrophysical, nuclear and particle Physics scenar- 
ios, the reached intrinsic radiopurity, the used approach, the specific performances and 
operating conditions, the large collected exposures of the former DAMA/Nal and of 
the present DAMA/LIBRA set-ups, have offered an unique possibility of an effective 
model independent investigation. 

Model dependent considerations will be presented in later publication specifically 
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devoted to this aspect. Just few arguments for some illustrative purposes are given in 
Appendix A. 

The collection of a larger exposure with DAMA/LIBRA (and with the possible 
DAMA/lton, which is at R&D stage) will also allow the improvement of corollary 
information which can be derived on the nature of the candidate particle (s) and on 
the various related astrophysical, nuclear and particle Physics scenarios, and the inves- 
tigation with very high sensitivity of the other DM features and second order effects 
as well as of several rare processes other than DM. 



6 Appendix A 

As in the past (see e.g. [4, 5, 6, 7, 8, 9, 10, 11]), corollary investigations can also 
be pursued - on the basis of the cumulative 8.2 a C.L. model-independent result by 
DAMA/Nal and DAMA/LIBRA - on the nature of the DM candidate particle and 
on related astrophysical, nuclear and particle Physics scenarios. As widely discussed 
elsewhere, these investigations are instead model-dependent and - considering the large 
uncertainties which exist on the astrophysical, nuclear and particle physics assumptions 
and on the theoretical and experimental parameters needed in the calculations have 
no general meaning (as it is also the case of exclusion plots and of the DM particle 
parameters evaluated in indirect detection experiments). 

Complete model dependent analyses, to update the allowed regions in various sce- 
narios and to enlarge the investigations to other ones, will be presented elsewhere. 
Here, we just remind that many astrophysical, nuclear and particle physics scenar- 
ios and many DM particle candidates exist. Just to offer some naive feeling on the 
complexity of the argument, we show in Fig. 21 the experimental S m values of Fig. 
9 with superimposed the expected behaviours for some DM candidates in few of the 
many possible scenarios and parameters values (see Table 6). In fact, despite the 
behaviour of the S m values can be effective a posteriori for template purpose, very 
accurate results on corollary model dependent quests at given C. L. should be evalu- 
ated by applying the maximum likelihood analysis in time and energy of all the events 
(as described elsewhere), which offers efficient and complete data analyses accounting 
for all the experimental information carried out by the data and, when of interest, for 
priors. 

It is worth noting that an increase of the exposure and a possible lowering of 
the used 2 keV threshold will improve the discrimination capability among different 
astrophysical, nuclear and particle Physics scenarios. 

Let us, finally, note that results obtained with different target materials and/or 
different approaches cannot intrinsically be directly compared among them even when 
considering the same kind of candidate and of coupling, although apparently all the 
presentations generally refer to cross section on nucleon. 

For completeness, we also further note that no experiment exists whose result can 
be directly compared in a model independent way with the ones by DAMA/Nal and 
DAMA/LIBRA. Thus claims for contradictions are arbitrary, in fact, e.g.: 1) the 
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Figure 21: Few examples of expected behaviours for few of the many possible scenarios 
(see Table 6), superimposed to the measured S m ,k values of Fig. 9. The shown 
behaviours have not been obtained by the maximum likelihood method (see in our 
quoted literature) and are shown just for illustrative purpose; they all give practically 
about the same C.L.. As mentioned in the text, the full treatment of the data by 
maximum likelihood method to update the volumes/regions allowed at given C.L. by 
the cumulative DAMA/Nal and DAMA/LIBRA data for the considered scenarios will 
be presented elsewhere following the full analysis method of refs. [4, 5, 6, 7, 8, 9, 10, 11]. 



others arc insensitive to the annual modulation signature; 2) the others use different 
target materials; 3) the plots they show and those they attribute to others are built 
at a "single cooking" without accounting at all for the existing experimental, theo- 
retical and phenomenological uncertainties and for the existing alternative choices; 4) 
DAMA/Nal and DAMA/LIBRA have a favoured sensitivity with respect to others in 
several scenarios; moreover, scenarios exist (see literature) to which the others are not 
only disfavoured with respect to the DAMA experiments, but even blind. Further- 
more, additional realistic limitations in those claimed model dependent sensitivities 
(just for "nuclear recoils" and a single assumed scenario and parameters set) arise so 
far e.g. from: 1) the unproved physical threshold with suitable keV source calibra- 
tions; 2) the energy scale extrapolated from higher energy; 3) unproved stability of 
the running parameters and of all the used "rejection" windows over long term at the 
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Table 6: Scenarios and parameters values used in Fig. 21; they have been chosen 
among the many possible ones [4, 5, 6, 7, 8, 9, 10, 11]. In the fourth column the 
considered Set - as in ref. [4] - of nuclear form factors and/or of nuclear quenching 
factors is reported. Here: i) asi is the spin independent point-like cross section; ii) 
<tsd is the spin dependent point-like cross section; iii) 9 is an angle defined in the 
[0,7r) interval, whose tangent is the ratio between the WIMP-neutron and the WIMP- 
proton effective SD coupling strengths, respectively [4]; iv) tuh is the mass of the 
LDM particle; v) A is the mass splitting [11]; vi) g aee is the bosonic axion-like particle 
coupling to electrons. For the cross sections of the LDM particle sec rcf. [11] and for 
the channeling effect see ref. [9]. 
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needed precision; 4) insensitivity to candidates giving part (even WIMPs) or all the 
signals in electromagnetic form; 5) marginal exposures; 6) unproved determination of 
the efficiencies in each one of the many data handlings, they apply at the needed 
level of claimed precision (a control of systematics at level of 10~ 4 - 10~ 8 is required); 
7) disuniformity in the detector response e.g. in two-phase liquid Xenon detectors; 
etc. Moreover, they generally quoted in an uncorrect, partial and unupdatcd way the 
implications of the DAMA/Nal model independent result. Some arguments have been 
addressed e.g. in rcf. [4, 5, 6, 7, 8, 9, 10, 11, 35] and in some literature. On the other 
hand, whenever there might be in future some correct claim for exclusion in one or 
more particular astrophysical, nuclear and particle physics model framcwork(s) and as- 
sumed parameters set (but correctly accounting both for experimental and theoretical 
uncertainties and for the implications of the DAMA model independent results in the 
considered scenario and assumptions), there will be still many other scenarios, param- 
eters sets and DM candidates which can explain the DAMA/Nal and DAMA/LIBRA 
model independent results to which other target materials and approaches are dis- 
favoured or even blind; thus, this never will exclude the model independent results of 
the DAMA experiments. In addition, whenever an experiment using the same identical 
target material and methodological approach would be available, as usual in whatever 
field of Physics a serious comparison would require - in every case - e.g. a deep in- 
vestigation of the radiopurity of all the part of the different set-ups, of their specific 
performances in all the aspects, of the detailed procedures used by each one, etc. 
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